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We examine the effect of an impurity on the nearby tunneling conductance in 

o 

an anisotropically-gapped superconductor. The variation of the conductance has 
pronounced spatial dependence which depends strongly on the Fermi surface location 
of gap extrema. In particular, different gap symmetries produce profoundly different 
spatial features in the conductance. These effects may be detectable with an STM 
measurement on the surface of a high-temperature superconductor. 
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Any well-developed theory of high-temperature superconductivity must predict the sym- 
metry of the energy gap. Some proposals, including the strongly-coupled phonon-mediated 
pairing which causes superconductivity in ordinary metals [[[]], yield isotropic or nearly 
isotropic gaps. The dominant characteristic of the pairing interaction in this model is retar- 
dation. In the low-temperature superconductors this is reflected in the frequency- dependence 
of the superconducting gap and clearly shows that phonons mediate the pairing. 

For a nearly half-filled Hubbard model, approximations involving the exchange of para- 
antiferromagnetic spin fluctuations give rise to an effective electron-electron interaction 
which has a characteristic momentum dependence, peaking near the antiferromagnetic wave 
vector 0H . A similar momentum dependence has been found in Monte Carlo simulations 
0. Such an interaction favors d x 2_ y 2 pairing and gives rise to a momentum-dependent gap 
which has four nodes. 

Gap-measurement techniques that were introduced and developed for use on electron- 
phonon superconductors measure accurately the frequency structure of the gap. Among 
these are the voltage dependence of the tunneling I(V) characteristic into the homoge- 
neous superconductor |J and the frequency-dependence of electromagnetic absorption |3|. 
If, however, the gap has strong momentum-dependence these measurements depend on some 
momentum- averaged value of the gap. Most of these probes |7||§ now indicate pronounced 
gap anisotropy or gapless superconductivity in the high-temperature superconductors. 

In the cuprate-oxide superconductors it would be useful to have experimental information 
on the momentum-dependence of the gap. Certain techniques exist which are designed to 
detect gap nodes on the Fermi surface. These include measurements of the low-temperature 
dependence of thermodynamic and transport properties, such as the specific heat || and 
the magnetic penetration depth ||10| . Here power-law dependencies on temperature typi- 
cally imply nodes, in contrast to the exponential behavior associated with a fully gapped 
Fermi surface. Measurements of these quantities proved extremely illuminating in studies of 
heavy- fermion materials Recent results |IT]] on YBa2Cii307 indicate a low-temperature 
penetration depth proportional to T, suggesting the existence of nodes. 
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At this time, however, only a few experiments exist that directly measure the g ap as a 
function of momentum on the Fermi surface. Angle-resolved photoemission (ARPES) on 
Bi 2 Sr 2 CaCu208 measures the magnitude of the gap [|12j with an angular resolution of about 
6° and an energy resolution of about 10 meV. A gap has been detected with ARPES only 
in Bi 2 Sr 2 CaCu 2 0g. Measurements of YBa 2 Cu307 — Pb SQUIDs may eventually determine 
the symmetry of the superconducting gap fl3fl , but do not provide detailed information on 



the momentum- dependence of the gap. In theory, measurements of linewidth changes of 



phonons with finite momentum have excellent angular and energy resolution ||14|| , and can 



measure the relative phase of the gap |LJ|. Such measurements have yet to be done. 

Here we propose that a scanning tunneling microscope (STM) study of the spatial vari- 
ation of the tunneling conductance |16j around an impurity can be used to probe the 



momentum-dependence of the superconducting gap A k . We will discuss in this Letter those 
features of the gap which produce the most profound changes in the local tunneling con- 
ductance: gap minima and maxima. A gap minimum or maximum on the Fermi surface 
produces strikingly large conductance variations along the spatial directions away from the 
impurity which are perpendicular to the Fermi surface tangent at the gap extremum. We 
numerically calculate the spatial structure for superconductors with (a) an isotropic gap, (b) 
a d x 2_ y 2 gap, and (c) an anisotropic but nodeless gap. The results are explained by analytic 
forms in various limits. The energies of gap maxima and minima should be discernible with 
a resolution better than 1 meV. 

The feasibility of observing the local tunneling density of states with atomic resolution 
has been established by STM tunneling conductance measurements near impurities and 
step edges on the surface of Cu(lll) [fT7|| . In addition, tunneling conductances have been 
obtained with near-atomic resolution on the surface of the high-temperature superconductor 
Bi 2 Sr 2 CaCu 2 8 @. 

We assume the impurity has three possible effects on the local environment in a supercon- 
ductor: an impurity potential 5e, a change in the superconducting gap 5A^ (which depends 
on the momentum k), and a change in the local mass Sm of the carriers. In this Letter each 



of these is modeled as a delta-function perturbation at an impurity site 0, perturbing the 
homogeneous Hamiltonian of the bulk superconductor. An impurity would actually have a 
finite-range effect on 8e, 5A\ l , and Sm, which could be as large as a superconducting coher- 
ence length. Here we will be concerned with conductance changes which persist at greater 
distances. 

Each of the perturbations contribute to the others; for example, introducing an impurity 
potential will have an effect of the same order on the gap magnitude. This in turn affects 
the conductance. The precise relationship between the parameters depends upon details of 
the microscopic mechanism. We will use Se and SA^ as phenomenological parameters in our 
model, and assume that they satisfy the self-consistent equations for the superconductor. 
The effect of a 8m can be absorbed into Se. We thus avoid the difficulty which results from 
attempting to start from "bare" impurity parameters. 

The effect of the impurity perturbations on the conductance under the STM tip can 
be calculated from the Green's functions of the homogeneous superconductor. These are 
written compactly in the Nambu formalism [113]: 



/G (x,x / ,u;) F (x,x',u;) \ 
G Q (x,x»= , (1) 

\F (x,x',u;) -G (x,x', -co) J 

where 

V V w - 4 - A k + «7 

*UX,X» = ^E^ (X ' X 2 2 Ak A 2^- • (2) 

Here is the single-particle energy in the normal material, Ak is the superconducting 
gap, x and x' are positions and V is the volume. The conductance at a position r and 
voltage V is the imaginary part of the upper- left entry in the matrix (|]): dl(r, V)/dV = 
— AImG (r, r, u = V). Here A consists of voltage and position-independent matrix elements 
and numerical factors. 

The Green's function with an impurity at to linear order in the perturbations is 

G(r, r, u) = G (r, r, u) + G (r, 0, uo) [<5eT 3 G o (0, r, u) + tfAnG^O, r, u)] , (3) 



where T\ and T3 are Pauli matrices, and 



^^-"'^•""^TO' < 4) 

Here 5A^ = 5Ad(k), where the maximum value of d(k) on the Fermi surface is 1. Since the 
Green's functions in Eq. (2) and Eq. (4) only depend on |x — x'|, the two position variables 
x and x' will be replaced by |x — x'| below. 

It is useful to express the conductance in terms of dimensionless quantities. We define 
Green's functions normalized by the density of states at the Fermi energy N*, such as 
9o — G a /N*. Then, in terms of the normalized couplings 5e = SeN* and 6 A = 5 AN*, the 
conductance at r is 



dl{r, V) 
dV 



Inu7 o (0, V)+83m {g 2 a (r, V) + / D 2 (r, V)} 



+<5AIm {^(r, V)f (r, V) + f' (r, V) 9o (r, V)} 



(5) 



where dl(r, V)/dV is normalized by the normal metal's conductance. 

We have plotted in Fig. 1 the normalized conductance (not including the bulk term) from 
Eq. (|]) for 6e — 1 and 6 A = for two gaps: (a) an isotropic gap, A(0) = A OJ (b) a d x 2_ y 2 
gap, A(0) = A o cos20. The bias V is set to 1.1 A Q . Strong angular-dependence of the signal 
in Fig. 1(b) is apparent. We have not plotted the conductance for the anisotropic s-wave gap 
A(0) = A o (0.55 + 0.45 cos 40), since at this bias voltage it appears essentially the same as 
Fig. 1(b). We will show later that it can be distinguished from the d x 2_ y 2 gap at lower bias 
voltages, near the anisotropic s-wave gap minimum. The Fermi surface is assumed circular, 
although the gap is considered pinned to the lattice. Our use of a cylindrical Fermi surface 
is not a significant assumption since, as will be explained below, the main contribution to 
any signal comes from small regions of the Fermi surface. 

In order to illustrate the strong angular-dependence of the conductance for the d x 2_ y 2 
gap, we show in Figs. 2(a) and (b) the conductance as a function of distance for two angles. 



The first angle, = 0, is the direction of a gap maximum. The second angle, = 7r/4, is 
the direction of a node. For comparison we show in Fig. 2(c) the conductance as a function 
of distance for an isotropic gap. 

We will now analyze the Green's function g Q to understand the pronounced spatial vari- 
ation in Figs. 1 and 2. It is straightforward to reduce the integrations for the Green's 
functions to angular integrals. Defining r = (x,y) = kp 1 p(cos (fi , sin0 o ), 
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Here Q(E) is the Heavyside step function. f (p,<f),uj) is given by a similar expression, with 
the parenthetical factors before the step functions replaced by A((j))/\u 2 — A 2 (0)|5. Similarly, 
g' and f' have the same form as g a and f a with an additional factor of d((f>) in the integrand. 

The largest contributions to the type of integral in Eq. (6) come from two occurrences in 
the integrand. The first occurrence is when the bias u is very close to the gap A(0). In this 
instance the denominator of the integrand becomes small. The second occurrence is when 
the curvature of the exponentiated factor becomes small. Since / d(f> exp(o;0 2 ) oc a -1 / 2 , when 
the curvature a is small the integral becomes large. To locate angles where the exponent has 
small curvature, we have expanded the ^-dependent factors of the exponent. The curvature 
is a minimum when A(0) is a minimum or maximum, and uu ~ A(0). The integration near 
these angles provides most of the variation in the tunneling conductance and we find the 
following form for the local contribution to the conductance in a d x 2_ y 2 superconductor: 

dI( P ,(f) ,V) -7T / V 2 



dV p cos (j) Q a \ V 2 — A 2 / 
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For brevity, the solution for — 7r/4 < O < 7r/4 is presented above. All other angles O can 
be mapped into this region. A<p is a p-independent angle and a is the curvature at the 
extrema, which is approximately given by 

1 2A 2 

a = °- r . 8 

2 v F k F {V 2 - A 2 )5 

The conductance in Eq. (0) has oscillations whose wavelength depends on the bias. This 
bias dependence should help screen out spurious effects. 

The isotropic-gap superconductor has no Fermi surface location where the curvature 
becomes smaller than elsewhere. Thus we find, as shown in Fig. 2, that the amplitude of the 
conductance variations is larger for the d x 2_ y 2 gap, in the spatial direction corresponding to a 
gap extremum, than the amplitude for the isotropic gap. The conductance for the isotropic- 
gap superconductor is given by the expression in Eq. (0) with the sin 2 (p sin O A0) / (p sin (fi ) 2 
and cos O factors absent and a = 1/2. 

Since the main contribution to the high-signal regions comes from near the gap extrema 
(where the curvature is least), it is reasonable to set uo = A Q = A(0) in these integration 
regions. Then f = 9o = j' = g' - The analytic work done on g Q above thus applies to 
all the Green's functions. We find numerically that plots of the change in the tunneling 
conductance due to a change in <5Ak are almost identical to those for a change in Si. 

In Fig. 3 we show the normalized conductance for the anisotropic s-wave gap and the 
d x 2 -y 2 g a P a t a bias V = .16A Q (this is just above the minimum gap for the anisotropic s-wave 
gap). The isotropic gap is not shown because there is no variation in its conductance due to 
the impurity at this bias. Figs. 3(a) and (b) show the change in the tunneling conductance 
for the anisotropic s-wave gap as a function of distance in the direction of the gap minimum 
and maximum respectively. The bias is close to the minimum, so there is a long-range 
enhanced signal only for (a). In contrast to (a) and (b), the change in conductance for the 
d x 2_ y 2 gap in the same directions, shown in (c) and (d), is small. 



Here we have shown that there is a characteristic structure in the local conductance 
at bias voltages near gap minima or maxima which can be used to probe the momentum- 
dependence of Ak- Thus one can obtain information about the symmetry of the gap. Beyond 
this, we believe that scanning tunneling microscopy of the local density of states around an 
impurity site can provide detailed information on both the momentum and frequency de- 
pendence of the gap. This would provide a means of determining both the momentum and 
frequency dependence of the interaction responsible for pairing in the cuprate superconduc- 
tors. 
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FIGURES 

FIG. 1. Normalized conductance for a unit perturbation to the site energy Se. The length 
scale is set by ftp 1 . Here we have taken a coherence length equal to lOfc^ 1 . The voltage is 1.1A , 
where A is the gap maximum. The range of the intensity is from -0.3 (black) to 0.3 (white), (a) 
isotropic gap, A(</>) = A Q (b) d x 2_ y 2 gap, A(</>) = A Q cos2</>. 

FIG. 2. Same function as Fig. 1. (a) Radial plot of the change in the conductance due to 
Se along the high-signal direction </> = for the d x 2_ y 2 gap. (b) Radial plot along the low-signal 
direction, <p = it/ 4. (c) Radial plot for the isotropic gap. The scale is the same for all three plots. 

FIG. 3. Normalized conductance for a unit perturbation to the site energy Se. Here 
the voltage, 0.16A o , is just above the gap minimum for the case of an anisotropic s-wave gap 
A(</>) = A o (0.55 + 0.45 cos 40). Radial plots of the change in conductance due to Se for the 
anisotropic s-wave gap are shown in the direction of the gap minimum (a) and gap maximum (b). 
Since the bias is close to the minimum, only (a) has a long-range enhanced signal. For comparison, 
the same functions for the d x 2_ y 2 gap are shown in (c) and (d). Compared to (a), all three of 
(b), (c) and (d) are small. Thus we find a strong dependence of the tunneling conductance on gap 
minima, similar to the dependence on gap maxima shown in Figs. 1 and 2. 
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